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SUMMARY 


To  develop  the  analysis  and  design  methods  for  robust  control  of  large  space 
structuie  sampled  data  stochastic  systems  ,  the  theory  of  multi-input,  multi-output  transfer 
function  matrix  with  the  z-transformation  is  used.  The  analysis  and  design  of  robust  control 
of  multivariable  discrete-time  feedback  systems  in  the  frequency  domain  are  studied.  The 
singular  value  theory  is  used  to  establish  robustness  stability  criterian  for  discrete-time  sys¬ 
tems  in  the  presence  of  parameter  urcertainties  or  unmodelled  dynamics  in  the  frequency 
domain.  The  robustness  stability  conditions  for  discrete-time  systems  with  additive  alteration, 
and  with  multiplicative  alteration  are  developed.  The  linear-quadratic-regulator(LQR)  and 
linear-quadratic-Gaussian(LQG)  theory  are  used  for  studying  the  loop  transfer  recov- 
ery(LTR)  of  discrete-time  systems,  and  the  robust  LQG/LTR  method  has  been  extended 
from  the  continuous-time  systems  to  the  discrete-time  systems.  It  has  been  proven  that  the 
LQG/LTR  method  is  also  valid  for  the  LQG  control  of  the  discrete-time  systems  with  the 
filtering  observer. 

As  an  application  of  the  LQG/LTR  technique  for  discrete-time  systems,  the  design 
of  reduced  order  optimal  digital  LQG  controllers  for  the  orbiting  flexible  shallow  spherical 
shell  system  is  considered..  Simulations  have  certified  the  12-dim.  reduced  order  controllers 
will  be  sufficient  for  the  optimal  LQG  control  of  the  orbiting  shallow  spherical  shell  system 
in  the  presence  of  unmodelled  dynamics.  The  performance  of  the  8-dim.  reduced  order 
LQG  controllers  for  the  shell  system  is  unacceptable.  The  6-dim.  reduced  order  controllers 
for  the  shell  system  will  result  in  the  severe  divergence  of  the  transient  responses. 

The  comparisons  between  the  digital  optimal  LQG  controller  with  the  filtering 
observer  and  predicting  observer  for  the  orbiting  flexible  shallow  spherical  shell  system  have 
been  made.  Based  on  a  comparison  of  the  robustness  recovery  properties,  the  system  per¬ 
formance  of  the  robust  control  system  with  the  filtering  observer  has  better  transient  response 
characteristics  than  the  LQG  robustness  control  system  with  the  prediaing  observer. 
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Introduction 


The  purpose  of  this  research  is  to  study  and  develop  the  analysis  and  design  meth¬ 
ods  for  robust  control  of  the  large  space  structure  sampled  data  stochastic  system  with  a 
specific  application  to  the  orbiting  flexible  shallow  spherical  shell  system. 

It  is  well  known  that  one  of  the  most  important  breakthroughs  in  multi-input, 
multi-out  feedback  system  theory  for  the  last  decade  is  the  deveU'pment  of  the  loop  transfer 
recovery  methodology  for  the  linear  quadratic  Guanssian  problem  which  is  called  LQG/ 
LTR.  Unfortunately,  the  previous  research  works  in  this  field  are  almost  all  for  the  continu¬ 
ous-time  system.  Therefore,  we  must  study  and  address  the  following  problems  before  ap¬ 
plying  the  LQG/LTR  technique  to  design  discrete-time  system  robust  controllers  for  the 
orbiting  large  space  structural  system. 

(1)  How  to  get  the  frequency  response  of  the  transfer  function  matrix  for  the 

discrete-time  systems  supposing  the  discrete-time  system  model  is  given 
in  the  time  domain  ? 

(2)  What  is  the  robustness  stability  condition  for  the  discrete-time  system  in 
the  frequency  domain  ? 

(3)  Is  there  loop  transfer  recovery  for  the  discrete-time  systems  7  How  to  prove 
it  ? 

!n  this  report,  the  first  three  Chapters  will  address  the  three  problems.  As  an 
application  of  the  LQG/LTR  technique  for  the  discrete-time  system,  Chapter  4  studies  the 
design  problem  of  the  reduced  order  optimal  digital  LQG  controller  for  the  orbiting  flexible 
shallow  spherical  shell  system. 

Chapter  6  considers  the  differences  between  the  digital  optimal  LQG  controller 
with  filtering  observer  and  predicting  observer  for  the  orbiting  flexible  shallow  shell  systen. 
The  general  conclusions  and  suggested  future  direction  will  be  given  in  Chapter  7. 


I 


I  The  Analysis  Method  in  the  Frequency  Domain  for  Discrete-time  Systems 

As  we  know,  a  lot  of  research  about  robust  control  has  been  completed  for  continuous¬ 
time  systems  and  the  Laplace  tran'Tormation  method  is  its  main  analysis  tool  in  the  frequency 
domain.  What  we  want  to  do  is  to  find  the  relationships  between  the  transfer  functions 
of  the  continuous-lime  system  and  the  di.screle-timc  system  in  the  frequency  domain  ,  so 
that  the  research  results  of  the  robustness  problem  for  the  continous-time  system  may  be 
applied  to  the  discrete-time  system. 

1  I  Discrete-time  function  and  equivalent  continous-time  ''unction 

It  is  supposed  that  g(t)  is  a  continous-time  function.  g(kT)  is  the  sampling  function 
ol  g(t)  at  t=kT  points  (k=l,2,...),  i.e.  g(kT)  may  be  looked  at  as  the  output  of  an  ideal 
impulse  sampler  for  which  the  input  is  the  continous-time  function,  g(i). 
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Fig.  1-1  The  discrete-time  function  and  its  equivalent  continous-time  function 


The  output  of  the  ideal  impulse  sampler  can  be  written  mathematically  as  follows: 


OO 


g^(  t  )  =  I  g  (  t  )  8(  t-kT 
k 


OO 

=  Z  g  (kT)  5(  i-kT) 
k 


(1-1) 


Then  g^(t)  can  be  called  the  equivalent  continou''-time  function  associated  with  the  discrete- 
time  function  g(kT)  and  its  Laplace  transformation  i^ 
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oo 
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(1-2) 
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G"(  s  )  =  L{  g"(  I  )  }  =  J  Eg  (kT)  l-kT)  c  ''  cll  =  E  g  (kT)  e''*" 

k  k 

The  Z-iransformalion  of  g(k7)  is 

00 

G(  7.  )  =Z{  g(  kT  )}  =  E  g(  .^T  )  7-'^  (1-3) 

It  is  evident  that  the  relationship  between  tne  Laplace  transformation  of  the  equivalent  conti- 
nous-time  function  g*(t)  associated  with  the  discrete-lime  function  g(kT),  and  the  Z-trans- 
formalion  of  the  discrete-lime  function  g(kT)  is 

G(  z  )  =  G"(  s  )|  or  G^(  s  )  =  G(  z  )| 

s=(ln7)/  T  7=e®^  (1-4) 


It  will  be  proven  that  this  relationship  (1-4)  will  still  be  true  for  the  transfer  functions  be¬ 
tween  the  discrete-time  system  and  the  equivalent  continous-time  system. 

1.2  The  relationship  between  the  input  and  output  for  discrete-time  systems 
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Fig.  1-2  The  relati  mship  between  input  and 
output  fo  the  discrete-time  system 


It  is  supposed  that  the  ir  pul  and  output  for  the  discrete-lime  systems  is  u(kT),  and 
y(kT),  respectively.  The  relationship  between  the  input  and  output  in  the  time  domain  is 


n- 1  oo 

y(  nT)  =  E  g  ((n-k)T)  u(T  k  )  =  E  g  ((n-k)T)  u(T  k  )  (1-5) 

k  k 


(  because  g((n-k)T)=r  when  k  >  n-1  ) 
where  g  (nT  )  is  an  impulse  respon.se  function  of  the  discrete-time  system. 
The  Z-transformaiion  of  the  output  y(nT)  is 
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oo  oo  oo 

Y(z)  =  S  y(nT)  zr"  =  E  S  g  ((n-k)T)  u(kT) 
n  n  k 

oo  oo 

=  E  E  g  (mT)  u(kT)  7.-^'"^''^  (n-k=m) 

m  k 

oo  oo 

=  E  g(mT)  E  u(kT)  7.-''  (  g(mT)=0  ,  m<0  ) 

m  k 

=  G(z)  U(z)  (1-6) 

where 

oo 

G(z)=E  g(  mT)  7.-"’=  Z{  g(mT)} 

m 

oo 

U(z)  =  E  u(  kT  )  z-*'  =  Z{u(kT)} 
k 

1.3  The  relationship  between  the  input  and  output  or  the  equivaient  continuous-time 
system 


Fig.  1-3  The  relationship  between  input  and  output  for 
the  equivalent  continoas-time  sylcm 


It  is  assumed  that  the  input  and  output  for  the  equivalent  C(  niinious-time  system 
is  u^(t),  and  y^(t),  respectively. The  relationship  between  the  input  ard  output  in  the  time 
domain  is: 


where 


OO 

y^(  t  )  ^  /  gU  I  -  T  )  U*(  T  )  d  T 
0 


(1-7) 


OO 

u*(  I  )  =  u(  t  )  8r(  t  )  =  E  u(  kT  )  6(  t  -  kT) 

k 


OO 


g"(  t  )  =  g(  t  )  8r(  t  )  =  E  g(  kT  )  8(1  -  kT) 

k 
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00 

y^(  t  )  =  y(  t  )  5t(  I  )  =  S  y(  kT  )  6(  t  -  kT) 

k 

oo 

8r  (  I  )  =  S  8  (  I  -  kT  )  the  impulse  sampler 

k 

The  Laplace  transformation  of  (1-7)  can  be  expressed  as; 

OO  00  OO 

Y*(  s  )  =  J*  y^(  t  )  e  dt  =  J*  J*  g*(  t-T  )  u^(  t  )  e“®’  dt  dr 
0  0  0 
oo  00 

=  f  f  g^(  \  )  u^(  T  )  d\  dr  (  t-T=\) 

0  0 

oo  oo 

~  f  ^  f  ^  (  g*(k)=0  ,  k<  0) 

0  0 

=  G"(  s  )  U"(  s  )  (1-8) 

where 

OO 

G"(  s  )  =  /  g"(  \  )  e-*^  dv  =  L(g"(\)} 

0 

00 

u*(  S)  =  f  V*(\)  e-’^  d  .  =  L{u^(X)} 

0 

oo 

Y"(  s  )  =  /  y*(\)  e-*^  d\  =  L{y"(X)} 

0 

Because 

oo  oo  OO 

G^(  s  )  =  /  g  *(\)  e-*^  d\  =  /  E  g(  kT  )  8(  \-kT)  e’*^  d\ 

0  Ok 

oo  oo  oo 

=  E  /  g(  kT  )  e-*^  8(\-kT)  d\  =  E  g(kT)  6-“^* , 

k  0  k 

then  in  a  similar  way  one  can  obtain 
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OO  00 

U+(  s  )  =  /  u^(  t  )  e-*'  di  =  S  u(  kT)  e"'"' 

0  k 

OO  OO 

Y^(  s  )  =  f  y*(i)  e-*'  dt  =  Z  y(  kT)  6''''^' 

0  k 

Based  on  the  definition  of  the  Z-transformation  and  comparison  of  (1-6)  and  (1-8). 
we  have 

G"(  s  )  =G(z)  I  =  G(e*f) 
z=e*^ 

U^(  s  )  =U(z)  I  =  U(e*>)  (1-8) 

z=e®^ 

Y"(  s  )  =Y(z)|  =  Y(e^') 

z=e®^ 

Because 

^(s+jnW  )T  -^sT  gjnCn)  r=0sT 
(Os  =2Tr/T  T  —  sampling  period 

we  have 

G^(  s +jn(Os)  =  G"(  s  )  (1-9) 


Fig.  1-4  The  Frequency  Domain  for  the  Discrclc-time  System 
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The  following  transfer  function  properties  are  limited  to  a  particular  frequency  range,  indi 
caled  by  the  primary  strip  in  Fig.  1-4 

G"(  jw)  =  G(  )  (Ws  /  2  <  ca  <  cOs  /  2)  (I-IO) 

G^(jca)  =  G(ej‘’)  (-tt  <  0  <  tt)  (1-11) 

1.4  The  relalionship  between  the  input-output  for  the  equivalent  continuous-time 
system  and  the  discrete-time  system 

Because  the  output  of  the  equivalent  continous-time  systems  is 

OO 

y^(  l  )  =  L-'{  Y"(s)  }  =  L-'(G^(s)  U"(s)}  =  /  g*(  t-r)  u*(t)  dr 

0 

OO  OO  OO 

=  /  Z  g(kT)  &(  I-  T-kT)  Z  u(mT)  8(T-mT)  dr 

Ok  m 

OO  OO  OO 

=  Z  Z  g(kT)  u(mT)  /  8(  t-r-kT)  S(  r-mT)  dr 

km  0 

OO  OO 

=  Z  Z  g(kT)  u(mT)  8(  l-mT-kT) 

k  m 

OO  OO 

=  Z  Z  g(kT)  u(mT)  8t  t-nT)  (  m+k=n  ) 

k  m 

OO  OO 

=  Z  Z  g((n-m)T)  u(mTi  &(  t-nr)  (  g((n-m)T)=0,  n<m) 
m  n 

OO 

=  Z  y(nT)  8(  t-nT) 
n 

where 

OO 

y(nT)  =  Z  g  ((n-ni)T)  u(mT) 
m 

It  shows  that  th  :  relationship  between  the  output-input  of  the  equivalent  continuous-time 
system  is  simila  to  that  of  the  discrete-time  system. 


1.5  Summary 


(1)  If  the  transfer  function  of  the  discrete-time  system  is  given  by  Gfz),  (hen  (he 
transfer  function  of  its  equivalent  continuous-time  system.  G^(s),  can  be  expressed  a''  fol 
lows: 


G"(s)  =  G(z)| 
7,=e 


sT 


G’^Ow)  =  G(ei^)  (  -TT  <  9  <  7T  ) 

=  Gfe^^i^/^s)  (  -Ws  /2  <  u)<  a),/2  ) 

(2)  The  frequency  properties  for  the  discrete-time  system  can  be  described  in 
terms  of  frequncy  properties  for  the  equivalent  continuous-time  system. 
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2  Robustness  Criterion  in  the  Frequency  Domain  for  Discrete-Time  Systems 
with  Additive  Alterations  or  with  Multiplicative  Alterations 

2.1  The  basic  theorem  of  stability  for  IMIIMO  linear  feedback  system 


We  will  use  the  standard  nota  ion  of  input-output  stability  theory[7] 

a  =  some  Banach  space  of  function  x  :  T  X  with  ||  ■  || 

T  =  subset  of  the  real  numbers 
X  =  finite  dimensional  vector  space 
tte  =  {  X  :  PtX  e  a  for  all  t  gT  } 


(PTX)(t)  =  j 


*x(t) 


0 


t  ^  T 
t  >  T 


Jt.2'"  =  space  of  m-vcctor  function  on  T  with  integrable  Euclidean  norm 

5  :  a*  *  ae  =  identity  operator 
g  :  ae  ♦  a*  is  causal  if  Pt^Pt  =  PtO  for  all  t  g  T 

A’  =  conjugate  transpose  of  a  complex  matrix  A 


II  Prgxi  -PTgx2  II 

II  5  II  =  Sup  _ 

A  Xi,X2  €  ae  ||PtXi-PtX2||  t  G  T 

PtXi  X  PtX2 

We  consider  the  feedback  system  depicted  in  Fig. 2-1.  Here  the  causal  operator  g: 
l2e'"'^X-2e'"  represents  the  plant  plus  he  any  compensator  that  is  used. 
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The  basic  feedback  equation  is 

(5  +  0)e  =  r  (2-1) 

and  the  basic  stability  question  is  whether  (  3  +  Q  exists,  is  causal,  and 

is  a  bounded  operator  in  the  sense  that  jj  (  5  +  Q  )”’||  <  oo  .  We  will  assume  that  the 
nominal  system  is  stable.  We  are  interested  in  whether  the  closed  loop  sysem  retains  these 
properties  when  subject  to  additive  (  Q  0  +  ^0  )  or  multiplicative  + 

perturbations  representing  uncertainty  in  the  dynamic  behavior  of  the  system.  The  following 
theorem  provide  the  basis  for  our  analysis  .[7] 

Theorem  1  :  .X  :  tte  tte  be  a  linear  causal  operator,  and  suppose  exists,  is  causal 
and  m~'ll  <  oo  •  Then  if  AU  :  is  a  causal  operator  satisfying  11au||  < 

A  A 

oo  and  if 

||.X-’Aj11|  <  1  (2-2) 

A 

it  follows  that  ( jt+AU  )"^  :  ^  exists,  is  causal  and  has 


IK  .X+A,X  )~'||  <  oo  (2-3) 

A 

Theorem  2  :  (Desoer  and  Vidyasager  1975  [4J)  Let  the  operator  Q  :  i.je'"  — ^  Jfae"’  for 
T  =  (0,  oo]  be  defined  by 


00 

(Ox)(t)  =  /  Q(  t  -  T  )x(t)  dr  (2-4) 

0 

where  the  elements  of  the  impulse  response  matrix  Q(  t )  are  assumed  absolutely  integrable 
on  T  .  Then 


where 


lion  =  iioii  =  cT. 


(Tmax  =  max  max  ai(G(ju))) 
(i)>0,  l^i^m 


(2-5) 


and  where  ai(G(ju))  denotes  the  ith  singular  value  of  the  transfer  function  matrix  conre- 
sponding  to  Q  . 
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2.2  The  stability  conditions  Tor  the  MIMO  feedback  system  with  additive  alteration 
or  with  multiplicative  alteration 


Fig.  2-3  Feedback  System  with  Additive  Alteration 


Fig.  2-4  Feedback  System  with  Multiplicative  Alteration 


Where 


Qp  ;  lae" - ■t2e"'  the  causal  operator  of  the  plant 

Oc  :  Ize”'  — -  i-2e"’  the  causal  operator  of  the  compensator 

:  X-2e'"  ^  •t2e'"  the  causal  operator  of  the  additive  alteration 

AQm  :  i'2e'"  —*■  lit"'  the  causal  operator  of  the  multiplicative  alteration 
The  coresponding  transfer  functions  are  Gp.  Gc,  AGg,  AGm.  respectively. 

2.2.1  Feedback  system  with  additive  alteration 
For  the  feedback  system  with  additive  alteration.  Fig. 2-3.  we  have 
y=(0p  e  +  AQa  e)  =  (Qp  +  AQa)  e 
e=r-9c  y 

i.e.,  r=(  5  +  Oc  Op  +  Qc  Ag*  )e 

Applying  the  Theorem  1,  if 

IK  j  +  OcOp  r‘gcAg,ii  <  1  (2-6) 

A 

then 

(  i  +  9cSp  +  OcAOa  )"*  exist  and  is  bounded 

Applying  Theorem  2  to  (2-6),  and  considering  the  relationship  between  the  discrete¬ 
time  system  and  the  equivalent  continuous-time  system,  we  have  the  following  sufficient 
condition  for  system  stability; 

IK  i+Gc(e2'^i"^"»)Gp(e2^'“^"*))"'Gc(e2"'j"^"^)AGa(e^^j"'‘^')ll2  <  I 

V  CO  G  [  0,  10,72  )  (2-7) 

Because 

ct(A)  =  II  A  ||2  ,  ji(A)  =  (  II  A“*  ||2  )■' 
we  can  obtain  (2-8)  from  (2-7) 

II  (I+GcGp)“*GcAGa  II2  <  ||  (I+GcGp)"'Gc  ||2  ||  AGa  II2  <  1  (2-8) 
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If  Gc  '  exists,  then 


II  ^Ga  II2  <  ||(I+GcGp)"'Gc  II2  ^  =  ||(Gc  *  +  Gp)  '1I2  ' 
i.e., 

a[AGa(e2^j“^“*)]  <  2[Gc“  +  Gp(e2^j“^“*)] 

Vo)  E  [0.  a)s/2)  (2-9) 

or  we  have  from  (2-8) 

||(I+GcGp)-i||2  IlGclb  IIAGalU  <  1 
i.e., 

2(1  -Gc(e2’'j“'“»)Gp(e2^‘“'“*)] 

CT[AGa(e2^j“^“*)l  <  - 

2{Gc(e2’'j“^“*)) 

V  0)  E  [  0,  (Oj)  (2—10) 

2.2.2  Feedback  system  with  multiplicative  alteration 
From  Fig.  (2-4),  we  have 

y  =  gp(  5  +  AOm  )e 
e  =  r-Ocy 
i.e., 

r  =  (  3  +  Oc  Qp  +  9c  Qp  A9r 
Applying  Theorem  1,  if 

II  (  3  +  9cQp)  *  9c  9p  A^rn  II  f  (2~1  1) 

A 

then 


(  3  +  0c9p  +  9c9pA5m  )”‘  exists  and  is  bounded. 

Applying  Theorem  2  to  (2-11)  and  considering  the  relationship  between  the  discrete-time 
system  and  the  equivalent  continuous-time  system,  we  have  the  following  sufficient  condi¬ 
tion  for  system  stability: 
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II  (I+Gc-^Gw)Gpnja)))-'Gc"(ja))Gpnjw)AG,„njw)||2  <  l 


(OS 


=  27T/T 


II  (I+Gc(e2'"j“/“»)G(e2'"j“'“*))''Gc(e2^j“^“*)Gp(e2^j“^“*)ACm(e^^j“^“*)||2  <  1 

V  (d  6  {0.a)s/2) 

We  may  obtain  the  following  relationship  from  (2-  12) 

II  (I+GcGp)''GpGeAGm)||2  <  ll(I+GcGp)"’GcGp||2  llAGmll2  <  1 
i.e., 

afAGm]  <  l/a((I+GcGp)"'GcGp] 

If  (GcGp)'*  exists,  then 

CT[AGn,]  <  2[  (GcGp)-‘  +  I  1 
i.e., 

CTfAGmCe^'^'^'"*)]  <  a((Gc(e2^i“^“»)Gp(e2^j“^‘’ ^))-'  +  I  ] 

V  u  6  [0,  (0,/2) 


Because 


(I+A)+(H-A"*)  =  (I+A)(H-A“>) 


then 


(I+A)-'+(I+A-i)-‘  =  I 

$ 

||(I+A-»)-‘||2  -  II(I+A)-‘||2  ^  1 


i.e., 


l/ff[I+A-*]  ^  1  +  l/g[I+A) 


i.e.. 


2[I+A‘‘]  ^  g[I+Al/(l+s(I+A]) 


(2-12) 


(2-13) 


(2-14) 


-15) 
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Considering  (2-14)  and  (2-15),  we  obtain  (2-16)  if  (GcGp)  '  does  not  exist 


a[  I+Gc(e2’'i"^“*)Gp(e2^j“^“')] 

a[AGn,(e2^i“^“*)l  <  - 

I  +  g(I+Gc(e2^j“^"*)Gp(e2^j“'"*)l 

<  2((GcGp)"*+ I)  Vue[0,  a)s/2) 
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3  The  Linear-Quadratic-Gausslan(LQG)/Loop-Transrer-Recovery(LTR) 
Methodology  for  the  Discrete-time  System 

3.1  Introduction 

One  of  the  most  important  breakthroughs  in  muhi-inpul,  multi-output  feedback 
system  theory  for  the  last  decade  is  the  development  of  the  loop  transfer  recovery 
methodology  for  the  linear  quadratic  Gaussian  problem,  which  is  called  LQG/LTR.  The 
people  who  developed  this  techniques  at  first  arc  Kwakernaak. Doyle  and  Stein  [1,2| 

These  are  design  techniques  which  allow  tlie  excellent  robustness  and  sensitivity 
properties  of  optimal  state  feedback  schemes  to  be  almost  recovered  by  output  feedback 
schemes.  Although  this  was  the  original  motivation  for  the  development  of  these  tech¬ 
niques,  a  w'ider  and  more  important  aspect  of  them  is  that  they  simplify  the  use  of 
LQG  methodology,  allowing  practical  feedback  d  signs  lo  be  attained  with  a  reasonable 
amount  of  effort. 

Whereas  output  feedback  design  via  LQG  methods  usually  requires  the  specification 
of  two  pairs  of  matrices,  namely  a  pair  of  cost-wcighiitig  matrices  and  a  pair  of  noise 
covariance  matrices,  the  asymptotic  recovery  approach  -equires  only  one  of  these  pairs 
to  be  designed,  w'ith  the  other  pair  being  assigned  values  according  to  an  automatic 
procedure.  This  results  in  a  tremendous  reduction  in  the  ( omplexity  of  the  design  process. 
Consequently,  it  is  of  great  importance  to  obtain  an  analogous  procedure  for  discrete¬ 
time  systems.  Unfortunately,  the  previous  research  work  in  this  field  is  almost  all  for  the 
continuous-time  system  (except  Maciejowski  who  is  with  Cambridge  University,  England' 
The  difficulty  of  the  robustness  recovery  problem  (or  .sensitivity  recovery)  for  the  discrete¬ 
time  system  is  that  the  filter  gain  (or  control  gain)  for  tlie  discrete-lime  sy.stem  is  finite, 
but  the  filter  gain  (or  control  gain)  for  the  continuous-lime  system  is  infinite,  so  robustness 
(or  sensitivity)  recovery  can  not  be  obtained  by  simple  reduction  similar  to  the  case  for 
the  continuous  lime  svstem. 

Maciejowski's  work(10|is  only  about  the  development  of  the  sensitivity  recovery  of 
LQG  control  wi  h  the  filler  observer  for  the  discrcic-lime  system.  But  with  discrete-time 
sy.stems,  the  duality  with  the  optimal  control  is  only  for  the  predicting  observer, not  for 
the  filter  observer.  Therefore,  the  robustness  recoverv  for  the  discrete-time  system  can 
not  be  obtained  by  simply  using  the  duality  principle. 

Our  research  work  is  to  develop  techniques  for  ensuring  robustness  recovery  for  discrete- 
time  systems. The  specific  developments  are  as  follows. 


Fig. 3-1  Block  Diagram  of  ihe  LQG  Digital  Control  with  Filtering  Observer 

where  the  transfer  functions  of  the  plant  and  componsenior  are  Gp(7.)  and  Gc(7.),  they  are 
as  follows; 

Gp(7)=C(7l-A)-’B 
Gc(7.)=zKc(zI-(l-KfC)  (A-BKc)  )■’ K, 
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Fig. 3-2  Block  Diagram  of  the  LQG  Digital  Control  with  Predicting  Observer 

where  the  compensator  transfer  function  for  the  predicting  observer,  Gcp(7.),  is  as  follows 
as: 

Gcp(z)  =  Kc(zI-Ak)-'Kp 
Ak  =  A  -  BKc  -KpC 
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3.2  LQG/LTR  Methodology  for  Discrete-time  Systems 
3.2.1  The  statement  of  the  LQG  problem 

U  is  given  that  the  system  state  equation  is 

x(k+ 1  )=Ax(k)+Bu(k)+L^(k)  (3-1) 

The  measurement  equation  of  the  system  is 

y(k)=Cx(k)+p,Iq(K)  (3-2) 

The  control  output  equation  is 

yc(k)=Hx(k)  (3-3) 

It  is  assumed  that  the  statistical  properties  are  given  by 
E{e(K)C(K)r)  =  I  E{q(K)Ti(K)*'}  =  I  (3-4) 

then 


E{(U(k))(U(k))T}  =  Ltr  =  Q  (3-5) 

E{(p,Iti(k))(p,Iti(k))T}  =  p,2t  =  r  (3-6) 

where 

L  noi.se  input  matrix 
p  parameter  of  measurement  noise 
H  control  output  matrix 

p  parameter  of  control  weighting  matrix 

The  objective  of  the  LQG  problem  is  to  find  a  controller  depending  only  on  y(k),  u(k) 
(k=l .2,3,...),  to  minimize  the  performance  index,  J,  where 

oo 

J  =  E{  Z  (xT(k)Qx(k)  +  u'(k)ftu(k)  }  (3-7) 

k=0 

where 


H^H  =  Q  p^I  =  R 


(3-8) 


It  is  well  known  that  if  the  sys  em  (A  B  H)  is  controllable  and  observable  and  the  system 
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(A  L  C)  is  controllable  and  observable(  the  condiii  ms  may  be  reduced  to  stabili/.able  and 
detectable  for  the  time-invariant  system),  then  the  tiosed-loop  system  of  the  LQG  optimal 
controller  is  asymptotically  stable.  It  should  be  pointed  out  that  the  free  parameters  of  this 
problem  are  L  (noise  input  matrix),  |jl  (intensity  of  the  observational  noise),  H  (cotnrol 
output  matrix),  and  p(control  weighting  factor  in  the  performance  index).  In  typical  LQG 
applications,  these  parameters  are  assigned  a  priori  physical  significance  (e.g.,  process  noi.se, 
sensor  noise,  controlled  variables,  and  control  weights).  Then  the  solutions  for  the  LQG 


problem  are  : 

u(k)  =  -Kc  x(k/k)  (  For  the  filter  observer)  (3-9) 

or 

u(k)  =  -Kc^(k/k-I)  (For  the  predicu  r  observer)  (3-10) 

where 

x(k/k)=)5(k/k-l)+K((y(k)-9(k/k-l))  (3-11) 

x(k+l/k)=Ax(k/k-l)+Bu(k)-»-Kp(y(k)-y(k/k-l))  (3-12) 

y(k/k-l)=C«(k/k-l)  (3-13) 

Kp=AK,  (3-14) 

K{=P,CT(CPeC'r+R)-'  (3-15) 

Pe=APeA'‘-APeCT(R+CPeCT)''CPeA^+Q  (3-  1 6) 

and 

Kc=(ft+B^PB)-‘BPA  (3-17) 

P=A^  PA-A'^PB(ft+B'rpB)-'B'fPA+C  (3-18) 


The  block  diagrams  of  the  optimal  LQG  digital  control  with  filter  observer  and  with  predict¬ 
ing  observer  are  shown  in  Fig.  3- land  Fig.  3-2,  respectively. 


3.2.2  The  development  of  the  transfer  function  matrices 

From  (3-11)  and  (3-13),  we  have 

x(k/k)=(l-K,C)x(k/k- 1  )+K,y(k)  (3-19) 

Considering  (3-9),  (3-12),  and  (3-13),  wt  have 
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x(k+  l/k)=(A-BKc)  (I-K,C)x(k/k- 1  )+(A-BKc)  Kfy  (k) 
Let  x-(7.)=Z{x(k/k-l)}.  x(z)=Z{x(k/k)},  y(7.)=Z{y(k)} 


(3-20) 


(3-21) 

The  Laplace  transformation  of  (3-20)  can  be  written  as  follows 

x-(7.)=fzI-(A-BKc)(I-K,C)l-‘(A-BKc)Kry(7)  (3-22) 

The  Laplace  transformation  of  (3-19)  is 
x(7)=(I-KfC)x“(z)+Kry(7,) 

=(I-Kf)  [zI-(A-BK.)  (I-KfC)  |-’  (A-BKc)  Kfy  (z)+K,y(z) 

=f(I-KfC)[zI-(A-BKc)(I-KfC)r'(A-BKc)+I}K,y(z)  (3-23) 

Because  of 

x(zl-yx)"'y+l=7,(zl-xy)"‘  (3-24) 

then  (3-23)  can  be  written  as  follows: 

x(z)=z(zI-(l-Kf(.:)(A-BKc)]-'K,y(z)  (3-25) 

Let  G,-(z)  be  the  transfer  function  of  the  plant; 

Gc(z)  be  the  transfer  unction  of  the  compensator  with  filtering  obserber; 

Gcp(z)  be  the  transfer  function  of  the  compensator  with  predicting  observer. 

Then  the  compensator  transfer  function  for  the  filtering  observer,  Gdz),  can  be  obtained 
from(3-9)  and  (3-25)  as  : 

Gc(z)=zKc(zI-(I-K,C)(A-BKc)r’K,  (3-26> 

The  Gp(z)  can  be  obtained  from  (3-T)and  (3-2): 

Gp(z)=C(zI-A)-'B  (3-27) 

The  compensator  transfer  function  for  the  predicting  observer,  Gcp(z),  can  be  obtained  from 
(3-l0),(3-12),and(3-l 3)  and  is  expressed  : 

Gcp(z)-Mzl-Ak)  'Kf  (3-28) 

where 

Ak=A-BKc-KpC  (3-29) 
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3.2.3  Robustness  recovery  and  sensitivity  recovery 

It  is  well  known  that  the  multivariable  linear-quadratic(LQ)  optimal  regulators  have 
impressive  robustn'iss  properties,  including  guaranteed  classical  gain  margins  of  -6  db  to 
+0O  db  and  phase  nargins  of  +60°  in  all  channels.  The  result  is  only  valid,  however,  for 
the  full-state  case  If  observer  or  Kalman  filters  are  used  in  'he  implementation,  no  guaran¬ 
teed  robustness  properties  hold.  The  robustness  recovery  means  that  if  the  measurement 
noise  parameter,  p.,  approaches  zero,  then  the  loop  transh  r  function  of  the  LQG  control 
at  the  input  loop-breaking  point  1,  will  approach  t  le  loot)  transfer  function  of  the  LQR 
control.  The  sensitivity  recovery  means  that  the  loop  transf:r  function  of  the  LQG  control 
at  the  output  loop-breaking  point  2,  will  approach  ih :  loop  ransfc'^  function  of  the  Kalman 
filter  when  the  weighting  parameter,  p,  approaches  zero [  -3].  ^v'e  will  prove  that  these 
results  are  also  true  for  discrete-time  systems. 

The  following  facts  can  be  found  from  Fig.  3-1: 

(a)  The  loop  transfer  function  obtained  by  breaking  the  LQG  loop  at  point,  1’.  is  the  LQR 

loop  transfer  function  Ti(z),i.e.,Kc(zI-A)"'B; 

(b)  The  loop  transfer  function  obtained  by  breaking  the  LQG  loop  at  point,  1 ,  is  GcGp: 

(c)  The  loop  transfer  function  obtained  by  breaking  the  LQG  loop  at  point,  2’,  is  the 

Kalman  filter  transfer  function,T3(z),  i.e.,  C(zI-A)”’Kp: 

(d)  The  loop  transfer  function  obtained  by  breaking  the  LQG  loop  at  point,  2,  is  GpGc 


If  the  discrete-time  system  also  has  the  properties  of  the  loop  transfer  recovery  similar 
to  the  case  for  the  continuous-time  system,  then  the  following  relationships  should  be  satis¬ 
fied. 


Limt  GcGp  =  Kc(zI-A)“'B  (Robustness  recover”) 

and  0 

Limt  GpGc  =  C(zI-A)'*Kp  (Sensitivity  recovery) 

p-  0 

Before  beginning  the  formal  proof,  a  lemma  will  be  used: 

Lemma  (Shaked[lll)  :  If  det(CB)  ^  0,  and  system  (3-l)-(3-3)  is  of  minimum  phase, 
then  the  Kalman  filter  gain  Kp  determined  by  (3-l4)-(3-16)  will  be 
AB(CB)“'  when  the  variance  parameter  of  the  observational  noise,  p,, 
approaches  to  zero,  i.e., 

Limt  Kp  =  AB(CB)-' 

0 


(3-30) 

(3-31) 
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This  is  the  simplest  case  of  Shakccl’s  much  more  genera!  result.  Because  Kp=AKf 
and  null(A)=0,  therefore  Limt  Kf-B(CBr' 

Theorem:  If  the  open  loop  transfei  function  of  the  system  (3-l)-(3-3),  has  no  (finite) 
zero  in  {  z  :  z  >  1  }  and  det(CB)  ^  0.  then 

Limt  GoGp  =  Kc(zI-A)"'B  (3-32) 

IX— 0 


Proof: 

Let  A(z)  =  GcGp  -  Kc(zI'A)-'B 
Limt  A(z)=Limt  GcGp  -  Kc(zf-A)''B 

=Limt  ^Kc[zI-(I-Kf)(A-BKc)r'KfC(zI-A)-'B-Kc(zI-A)-'B 

p->0 

=Kc{Limt  z(zI-(I-K(C)(A-BKc)l-'K,C-I}(Iz-Ar'B 

=Kc{z(zl-(l-S)  (A-BKc)  r'S-I}(zI-A)'’ B 
=Kc{z[zI-(I-S)Aj-‘S-I}fzI-A)-'B 
=Kc(zI-(I-S)A]-'{zS-[zHl-S)Aj}(zI-A)-'B 
=Kc[zI-(I-S  A|-'{z(S-I)+(I-S)A}(zI-A)-'B 
=-Kc  (zI-(  1  -S)  A  r '  (I-S)  (zI-A)  (zl-A)' '  B 
=-K,(zI-(F-S)Ar'(l-S)B 

=  0  (  Since  (I-S)B=0  )  (3-33) 

where 


S=r-B(C  B)-'C 

The  satisfaction  at  Eq.(3-32)  means  that  there  is  the  robustness  recovery  properties 
for  the  digital  LQG  control  with  filtering  observer,  hut  there  is  no  the  robustness  recoverv 
properties  for  the  digital  LQG  control  with  the  predicting  observer,  because  the  Gc  and  Gep 
is  different,  (3-32)  can  not  be  proven  by  using  Gep  instead  of  the  Gc 

In  similar  way  the  (3-31)  can  be  also  proven,  it  indicate  that  the  sensitivity  recovery 
is  also  true  for  the  digital  LQG  control  with  the  filtering  observer. 
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Note; 


1 .  The  condition,  det(CB)  t  0,  ensures  that  Gp(z)  has  the  maximum  possible  number 
(n-m)  of  finite  zeros,  and  the  minimum  possible  number  (m)  of  infinite  zeros. 
The  perfect  recovery  obtained  is  only  possilile  because  the  nonzero  poles  of  Gc 

cancel  the  (n-m)  finite  zeros  of  Gp(z),  and  the  m  origin  poles  of  Gp(z)  cancel 
the  m  origin  zeros  introduced  by  the  factor  z  in  (3-26). 

2.  The  mechanism  by  which  the  recovery  is  achieved  is  essentially  the  same  as  in 
the  continuous-time  system  case  :  the  compensator  cancels  the  plant  zeros  and 
possibly  sorre  of  the  stable  poles,  and  inserts  the  controller  (observer)  zeros. 
Clearly,  this  will  fail  if  the  plant  has  zero.s  outside  the  unit  circle,  since  the 
compensator,  Gc,  guarantees  internal  stability. 
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4  Design  of  Reduced  Order  Optimal  Digital  LQG  Controller  for  the  Orbiting 
Flexible  Shallow  Spherical  Shell  System 


4.1.  Introduction 

Future  proposed  space  miss  ons  would  involve  large  inherentlv  flexible  systems 
for  use  in  communications,  radiometry,  and  in  electronic  orbital  based  mail  systems. 
The  u.se  of  very  long  shallow  dish-type  structures  to  be  employed  as  rcceivers/reflectors 
for  these  missions  has  been  suggested.  In  order  to  satisfy  mission  requirements,  the  pro¬ 
posed  LQG  digital  optimal  control  )f  the  shape  and  orientation  for  an  orbiting  shallow 
spherical  shell  are  also  studied  [1  |. 

Since  the  mathematical  systeni  model  is  inherently  of  high  order,  a  practical  con¬ 
troller  has  to  be  based  on  a  reduced  order  design  model.  It  is  the  purpose  of  this  chapter 
to  study  the  design  of  the  reduced  order  optimal  digital  controller  for  the  orbiting  shallow 
spherical  shell  system. 

One  of  the  most  important  breakthroughs  in  multi-input,  multi-output  feedback 
system  theory  for  the  last  decade  is  the  development  of  the  loop  transfer  recovery  meth¬ 
odology  for  the  linear  quadratic  Gaussian  technique,  which  is  called  LQG/LTR  [6-8). 
These  modern  techniques  are  very  useful  in  treating  unmodeled  dynamics  and  stochastic 
uncertainties  such  as  disturbances  and  sensor  noises.  As  we  know,  these  methods  have 
been  developed  only  for  the  continuous-time  system.  However,  in  practice,  observation¬ 
al  data  used  to  verify  the  orientation  and  shape  of  large  space  flexible  systems  will, 
in  general,  be  collected  on  a  sampled  basis  (  discrete-time  data  system).  In  order  to 
meet  the  requirements  of  design  for  the  discrete-time  data  system,  the  loop  transfer 
recovery  problem  for  the  discrete-time  system  has  been  developed  in  our  current  re 
search  work;  It  has  been  proven  that  the  robustness  recovery  property  for  the  LQG 
problem  is  also  true  for  the  discrete-time  system  if  the  open  loop  transfer  function  of 
the  system  has  no  (finite)  zeros  outside  the  unit  circle  and  det(CB)=0.  where  the  C  is 
the  observation  matrix,  and  B  is  the  control  influence  matrix. 

The  main  purpose  of  this  chapter  is  not  to  show  how  to  develop  this  method  theo¬ 
retically  in  detail  from  the  continue  u.s-time  system  to  the  discrete-time  system)  this  will 
appear  in  another  paper),  but  to  apoly  these  results  to  the  design  of  robust  reduced-ord¬ 
er  LQG  controllers  for  large  space  structural  systems  with  sampled  input  data. 

4.2  Mathematical  Models 

The  mathematical  model  of  an  isotropic  shallow  flexible  spherical  shell  in  orbit 
was  developed  in  Refs  [  14-16).  The  resulting  linearized  equations  of  motion  for  the 
rigid  rotational  and  generic  elastic  modes  were  developed  as  ; 


v|;  -  -  (1+m)^  =  Cx/Jx'"’wc^ 

$  +  4m,34>  -  (1-^3)4/  =  C./J.('’’a)c2  (4-1 ) 

10 

e  -  3^t20  -  (2I/Jy(°>)  Z  =  Cy/Jy(0’a)c2 

n=l 

€i  +  -  3)ei  +  (2I,(‘VMil)0  =  3I,0>/Mil  +  Ei/Mjlcoc^  (i  =  1,2.3....  10) 


where 


T  =  (Oct,  6i  =qi(t)/l  (i=l,2.3,...) 

The  derivative  in  Eq.(4-l)  is  with  respect  to  t. 

ij/.<t),0  =  yaw.  roll,  and  pitch  angles,  respectively,  between  the  undeformed  axes 
of  the  shell  and  the  axis  of  the  orbiting  local  vertical  /local  horizontal  system. 
qi(t)  =  modal  amplitude  of  the  ith  generic  mode. 

(Oc  =  orbital  angular  rate,  constant  for  assun  cd  circular  orbit. 

1=  charateristic  length  (the  base  radius). 

Mi  =  the  ith  modal  mass. 

=  principal  moments  of  inertia  of  the  undeformed  shell. 

Cx.  Cy,  Cz  =  the  components  of  external  torques. 

Xc  =  coordinate  of  differential  area  on  the  surface  above  the  base  plane. 

)ti.  1^2,  =  (Jz<°'-Jy(°^)/Jx(°’.  (Jx<°>-J7.(°’)/Jy(°\  (Jy^^^-Jx^°>) respectively. 

1,(0  =  /xc<Dx(0dv 
V 

(I>^(")  =  transverse  component  of  the  nth  modal  shape  function. 
cOn  =  natural  frequency  of  nth  mode. 

On  =  U)n/Wc 

The  mode  shapes  of  the  transverse  vibratioi  s  of  a  shallow  spherical  shell  with 
a  completely  free  edge  are  given  (13): 

‘t>x(">  =  Api  (  (lP"VRDM,pj^)CpjCP  +  .Ipl^tp,.  +  r)pjlp(tipj,  ?)}cosp(P+Po)  (4-2) 


where 


p=the  number  of  nodal  diameters(meridians 
j=the  number  of  nodal  circles 

Apj,  Cpj,  Dpj  =  the  shape  function  coefficient'; 

tJipj  =  frequency  parameter 


D  =  bending  .stiffness  factor 
R  =  radius  of  curvature  of  the 
shell 

Jp  =  the  Bessel  functions  of  the 
first  kind 

Ip  =  modified  Bessel  furi^tions 
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Ihble  4-1  The  first  en  natural  frequencies  of  the  shell 


n 

p 

j 

u)  ^  ( rad./.sec) 

im 

2 

0 

1.027715063 

0 

1 

1.027781054 

1 

1 

1.028163456 

2 

I 

1.029176311 

0 

2 

1.029465999 

1 

2 

1.031988029 

2 

2 

1.036220825 

0 

3 

1.036928081 

1 

3 

1.0445%445 

HB 

2 

3 

1.05.5618524 

It  is  assumed  liiat  12  point  actuators  are  used  to  control  attitude  and  shape  of 
the  shell  system.  The  placement  of  the  12  actuators  is  determined  by  using  the  degree 
of  controllability  of  the  con  rol  system  f  13). The  location  of  actuators  on  the  shell  are 
shown  in  Fig. 4-1  and  the  ollowing  Table  2. 


Table  4-2  Locations  of  the  12  actuators  on  the  sht  How  shell 


Actuator 

No. 

I  cKations  of  actuator  (  ?,p  ) 
i  P 

Directions  of  the  jets 
f,”  f,"  t," 

1 

0.84 

0 

0 

2 

0.84 

7T 

1  0 

0 

3 

0.57 

tt/4 

1  0 

0 

4 

0.57 

5Tr/4 

1  0 

0 

5 

0.40 

3tt/4 

1  0 

0 

6 

0.40 

7tt/4 

1  0 

0 

7 

0  28 

tt/2 

1  0 

0 

8 

0.28 

3tt/2 

1  0 

0 

9 

1  00 

Tr/4 

1  -.simT/4 

COS7T/4 

10 

1.00 

5ti/4 

1  .simT/4 

-cost:/ 4 

11 

1.00 

3tt/4 

1  -.sin3TT/4 

co.s3tt/4 

12 

1.00 

7Tr/4 

-COS37t/4 

It  is  assumed  that  two  Earth  sensors  and  two  Sun  sensors  arc  used  to  measure 
the  attitude  angle  between  the  local  vertical/the  vector  of  the  Sun  direction  and  the  roll 
axis,  pitch  axis  of  the  shell,  respectively  It  is  also  assumed  that  8  displacement  sensors 
are  used  to  measure  shell’s  transverse  displacement  parallel  to  the  shell’s  yaw  axis;  the 
displacement  sensor  are  considered  to  becolocated  with  the  actuators. 
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Fig.  4-1  Orbiling  Shallow  Spherical  Shell  System 


The  slate  equal  ons  of  the  system  uc  as  follows; 

X  =  AcX  +  BcU 


(4-3) 


where 

Equations  (3)  i;an  be  discretized  as  follows: 

X(k+1)  =  Ax(k)  +  Bu(k)  (4-4) 

where 

A=exp(AcT)  ,  T=sampling  time 

T 

B=f  exp{Ac(T-t)}Bcdt 
0 

The  discretized  observation  equation  is  as  follows: 

y(k)  =  Cx(k)  (4-5) 

where 

x(k)  e  R".  u(k)  e  R'".  y  e  R^  ; 

A  €  R'”'",  B  G  R"’"",  C  G  R^”*" 

n=26,  m=12,  r=12  for  the  full  order  design  model.  It  is  assumed  that  the  reduced 

order  design  model  may  be  divided  into  4  cases.  They  are  shown  in  Table  4-3. 


Tbble  4-3  Full  order  and  reduced  order  design  mexJels 


The  order  of  design  nuxlel 

MtxJes  included  in  mtxlel 

n 

m 

ca.se  1 

full  order 

rigid  +  ID  modes 

26 

12 

12 

ca.se  2 

18-dim  reduced  order 

rigid  +  D  mcxies 

18 

12 

12 

case  3 

12-dim  reduced  order 

rigid  +  3  modes 

12 

12 

12 

case  4 

8-dim  reduced  order 

rigid  +  1  modes 

8 

6 

6 

ca.se  5 

6-dim  reduced  order 

rigid  only 

6 

4 

4 

4.3  The  Design  of  the  Reduced  Order  LQG  Controller  and  the  Loop  Transfpr 
Recovery 

4.3.1  Loop  transfer  recovery  for  the  discrete-time  system 

It  is  assumed  that  the  system  state  equation  ,  measurement  equation,  control  output 
equation  and  performance  index  for  the  LQG  problem  are  as  follov  s; 


x(k+l)  =  Ax(k)  +  Bu(k)  +  L^(k) 

(4-6) 

y(k)  =  Cx(k)  +  plTi(k) 

(4-7) 

yc(k)  =  Hx(k) 

(4-8) 

00 

J  =  E  S  {  x’^(k)  0  x(k)  +  u^(k)  ^  u(k)  } 

(4-9) 

k=0 


where 


E{(L^(k))(L^(k)'T}  =  LL^  =  Q 

E{^(k)er(k)}  =  I 

(4-10) 

E  {(M.h(k))(M.lTi(k))T}  =  =  R. 

E{Tl(k)Tl(k)l}=  I 

(4-11) 

0  = 

ft  =  pM 

(4-12) 

It  is  well  known  that  if  the  system  (A.B.H)  is  controllable  and  observable, 
and  the  system  (A,L,C)  is  controllable  and  observable  ( the  conditions  may  be  reduced 
to  .stabilizable  and  detectable  for  the  time-invariant  system),  then  the  closed-loop  sys¬ 
tem  of  the  LQG  optimal  controller  is  asymptotically  stable. 

The  LQG  control  law  for  the  filter  observer  is  : 


u(k)  =  -Kcx(k/k) 

(4-13) 

where 

5^(k/k)  =  5J(k/k-l)  +  Kf  (v(k)-9(k/k-l)) 

^(k+l/k)  =  Ax(k/k-l)  +  Bu(k)  +  Kp  (y(k)  -^(k/k-l)) 
y(k/k-l)  =  C^J(k/k-l) 

Kp  =  AK{  filter  gain  matrix 

Kf  =  PeC'^CCPeC^+R)-' 

Pe=APeA^  -  AP,C'(R+CPeC')-'CP,A’+Q 
and 

Kc  =  (^  +  B'PB)”'BPA  control  gain  matrix 


^0 


P  =  -  A'^PB(R+B^PB)-'B*^PA+  Q 

The  diagram  of  the  input  and  output  in  terms  of  the  Z  transformation  is  as 


follows: 


Fig4-2  Block  diagram  of  the  optimal  digital  control  with  filtering  observer 


where 


d)(z)  =  (zI-A)"* 

(4-14) 

Gc(z)  =  zKc[zI-(I-Kf)(A-BKc)r‘Kf 

(4-15) 

Gp(z)  =  C<I>B 

(4-16) 

It  is  well  known  that  the  multivariable  linear-quadratic  optimal  regulator (LQR) 
has  impressive  robustness  properties.  But  if  the  observer  or  Kalman  filter  is  used  in 
the  implementation,  the  robustness  properties  of  the  system  will  be  degraded.  The  robust¬ 
ness  recovery  means  that  if  the  measurement  noise  parameter,  p,,  approaches  zero,  then 
the  loop  transfer  function  of  the  LQG  control  at  the  input  loop-breaking  point,  1,  will 
approach  the  loop  transfer  function  of  the  LQR  control,  i.e.,  robustness  properties  for 
the  LQG  control  will  be  the  same  as  for  the  LQR  control  when  the  noise  parameter, 
M-,  approaches  zero. 

In  terms  of  mathematical  notations,  this  property  can  be  stated  as  follows: 

If  the  open-loop  transfer  function  of  system  (6)- (8)  has  no  (finite) zeros 

in  {  z:  |z|>l  }  and  det(CB)xO,  then 

lim  GcGp  =  Kc<I)B  (4-17) 

p,  -W  0. 

where  the  KcOB  is  just  the  loop  transfer  function  for  the  LQR  control. 

In  the  same  way,  the  .sensiiivity  recovery  can  be  also  proven,  i.  e.. 


(4-18) 


lim  GpGc  =  COKp 

p— 0 

where  the  COKp  is  just  the  loop  transfer  function  for  the  Kalman  filter. 

It  is  evident  that  the  robustness  recovery  properties  for  the  LQG  control  may 
be  used  in  the  design  of  the  reduced  order  control  ers,  so  that  the  controllers  designed 
based  on  the  reduced  order  design  models  will  1  ave  strong  robustness  properties. 


4.3.2  The  compromise  between  the  performance  and  the  robustness 

It  is  well  known  that  the  basic  requirement*^  of  a  feedback  system  are; 

(1)  Stability:  bounded  output  for  all  bound  d  disturbances  and  bounded 

refer  ence  input; 

(2)  Performance:  small  errors  in  the  presence  of  dist  jrbances  and  reference  input; 

(3)  Robustness;  stability  and  performance  m  rintai  led  in  the  presence  of  model 

uncertainties. 

In  fact,  the  requirements  (1)  and  (2)  are  in  conflict  with  the  requirement  of  (3). 
In  order  to  meet  the  requirement  of  (1)  and  (2),  we  should  keep  the  sensitivity  func¬ 
tion  [8) 


S(z)  =  (I  +  Gp(z)Gc(z))-'  (4-19) 

as  small  as  possible  for  all  frequencies  z=exp{jG.r},  (  -a)s/2  <  to  <(0s/2  ,  T=2'n-/t0s). 
If  the  requirement  of  (3)  has  to  be  met,  we  shouli  keep  the  complementary  sensitivity 
function  T(z), 


T(z)  =  Gp(z)Gc(z)(I-l-Gp(z)Gc(z))-‘  (4-20) 

as  small  as  possible  for  all  frequencies  z=exp{j(ol },  (  -(Os/2  to  <(i)s/2  ,  T=2tt/{Os  ) 
However,  since 


S(z)  -»-T(z)  =  I  (4-21) 

they  cannot  be  made  small  simultaneously.  Rather  we  must  trade  off  the  size  of  one 
function  against  the  size  of  the  other  in  accordance  with  the  relative  importance  of  distur¬ 
bance/command  power  and  model  uncertainty  at  each  frequency. 

Tradeoff  between  transfer  functions  can  be  formalized  by  posing  them  as  func¬ 
tion-space  optimization  problems[8).  The  developments  in  [8]  are  for  the  continuous¬ 
time  system;  parallel  to  discussions  in  [8],  the  following  developments  can  be  obtained 
for  the  discrete-time  systems. 

The  first  thing  needed  for  the  optimization  problem  is  to  chose  a  convenient  criteri¬ 
on  of  smallness.  Consider 


.^2 


a|M|2  =  <  Tr(MM"j 


This  shows  that  a  matrix  M  will  be  small  if  Tr(MM’'l  is  small.  Using  this  latter  measure 
for  the  two  matrices,  S(e^“^)  and  T(e*“^),  adding  weights  W(e'“^)  to  trade  one  against 
the  other,  and  integrating  over  the  frequency  range  results  in  the  following  optimization 
problem. 

Given  the  plant  Gp(z),  weight  W(z),  and  sensitivity  and  complementary  sensitivity 
functions  defined  by  (4-19)  and  (4-20),  respectively,  find  a  .stabilizing  compensator 
Gc(z)  to  minimize 


a)s/2 

j  =  f  {Tr(SWW»S'' 

-Ws/2 


0)5/2 

+  Tr(TT"l}du)  =  / 
-0)5/2 


Tr|MM"|do) 


(4-22) 


where 

M=M(ei“'')  =  (S(ej“'')W(c'“'r),  T(ei“T)l 


In  mathematical  terms,  this  represents  an  H^-optimization  problem!  17].  Using  this 
approach  similar  to  that  in  [8],  it  can  be  proven  that  the  LQG  problem  for  the  discrete¬ 
time  system  defined  by  (4-6)-(4~l2)  may  be  converted  into  th)  equivalent  H^-optimi- 
zation  problem,  i.e.,  the  performance  index,  (4-9),  can  be  con'  erted  into  the  following 
formulation: 


J|qg 


where 


0)5/2 

(T'2'rr)  j  TrIP(e'"’ )P"(d“’)  jdw 
-0)5/2 


(T=2Tr/a)5) 


(4-23) 


P(ejo)r) 


HcDL-Hf  )BGc(I+GpGc)’'CcDL 


\-pGc(I+GpGc)''C<I)L 


-|xH(^BGc(I+GpGc)'' 

-pp,Gc(I+GpGc)~' 


(4-24) 


Comparing  (4-22)  and  (4-23),  it  is  now  apparent  thai  the  only  remaining  step 
needed  to  solve  (4-., 2)  is  to  fine  free  parameters  for  (4-24)  such  that  P(z)  reduces 
'o  IVl(z).It  is  easy  to  verify  that  t  le  following  choices  provide  the  desired  result. 


Choose  L  and  ja.  such  that  (Sensitivity  Rect'very) 


C<I)L/ja  ^  W(z)  (4-25a) 

and  let 

H=C  and  p  —  0  (4-25b) 

then 


/(I+GpGc)-'W  -GpGc(I+GpGc)-' 

P(z) - -  fJl 

0  0 


SW  -T\ 

0  0  / 

Note  that  for  each  nonzero  value  of  p,  the  LQG  solution  for  these  choices  produces 
an  appropriate  stabilizing  controller  which  is  H^-optimal  for  (4-23>.  Moreover,  since 
(4-23)  converges  to  (4-22),  a  sequence  of  decreasing  p  -values  produces  a  sequence 
of  controllers  which  optimizes  (4-22)  in  the  limit. 

It  is  also  easy  to  verify  that  the  following  alternative  to  (4-25)  produces  another 
useful  transfer  function  tradeoff: 


Choose  H  and  p  such  that  (Robustness  Recovery) 


j  (4-250 


and  let 
Then 


H<I>B/p  =  W(z) 
L=B  and  p,  —  0 


(4-26a) 

(4-26b) 


These  choices  accomplish  an  H^-tradeoff  between  the  sensitivity  and  complementary 
sensitivity  functions  at  the  input  loop-breaking  point  1  of  Fig.4-1  instead  of  at  the 
output.  Both  choices  (4-25)  and  (4-26),  will  be  referred  to  as  LQG/LTR. 

Eq.(4-I7)  (or  (4-18))  shows  that  the  optimal  loop  transfer  function  matrix  of 
a  minimum-phase  H^-problem  (4-22)  corresponds  to  the  loop  transfer  function  of 
the  LQR(or  Kalman  filter) problem  with  its  loop  breaking  point  l(or  2).  Moreover. 
Eq.(4-17)  (or  (4-18))  shows  that  the  sequence  of  LQG  solutions  generated  by 
(4-26)(or  (4-25))  converges  to  this  function(  ’’recovers”  this  function)  as  the  design 
parameter  tJL(or  p)  becomes  small.  These  results  may  be  listed  as  follows: 


Table  4-4  Robustn  ss  (Sensitivity)  Recovery 


Robustness  Recovery 

Sensitivity  Recovery 

Let  L  =  B,  r  >  m 

Let  H  =  C,  r  m 

Suppose  C(zI-A)"'L  is  min.  phase 

Suppose  H(zI-A)~'B  is  min.  phase 

When  p,  — “  0 

When  p  -^0 

then  T,(z)  —  Kc(zI-A)-‘B 

then  T2(z)  C(zI-A)-'Kf 

Relationship  between  LQ  regulator 

Relationship  between  Kalman  filter 

parameter  and  sensitivity  weighting 

parameter  and  .sensitivity  weighting 

H(zI-A)-’B/p  =  W(z) 

C(z!-A)•'L/^JL  =  W(z) 

A,  B,  H,  p  :  LQ  -  regulator  parameters 

A.  L,  C,  |JL  :  Kalman  filter  parameters 

W(z)  :  Sensitivity  weighting 

W(z)  :  Sensitivity  weighting 

Where 


XeR"-  UeR'".  YeR^ 

T|  (7-)=Gc(z)Gp(z):  The  transfer  function  at  the  input  loop-breaking  point  1 
T2(z)=Gp(z)Gc(z):  The  transfer  function  at  the  output  loop-breaking  point  2 
Kc(zI-A)”'B  :  LQR  loop  transfer  function 
C(zI-A)"*Kf  :  Kalman  filtei  loop  transfer  function 

The  significance  of  these  result  is  that  we  can  design  LQG  loop  transfer  functions 
on  a  full-state  feedback  basis  anc:  then  approximate  them  adequately  with  a  recovery 
procedure.  For  the  point  1,  'he  fu  l-state  design  must  be  done  with  the  LQR  equations 
and  recovery  with  the  Kalman  filtrr,  while  for  point  2,  full-state  design  must  be  done 
with  Kalman  filter  equation>  and  ecovery  with  LQR  equations. 

The  above  properties  suggest  a  tw'o-.step  approach  to  the  H^-optimal  problem 
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design: 

Step  1:  Design  a  LQ-regulator  (or  Kalman  (iker)via  (4-26a)(or  (4-25a)).  with 
desirable  sensitivity,  complemeniai  v  sensitivity,  and  loop  transfer  func 
tion. 

Step  2:  Design  a  sequence  of  a  Kalman  filier(or  LQ-rcgulaiors),  via  (4-26b) 
(or(4-25b)),  to  approximate  the  f'inction  in  step  I  to  whatever 
robustness  (or  accuracy)  is  needed 

Both  of  these  steps  are  easy  design  tasks  The  lirst  is  easy  because  the  LQ-rcguIa 
tor(or  Kalman  filter),  sensitivity,  complementary  seasitivity.  and  loop  transfer  function 
are  explicity  related  to  the  chosen  weights  W(ei‘^’^),  and  the  second  is  easy  because  it 
involves  only  repeated  solutions  of  the  algebraic  R  ccati  equations.  This  follows  by 
inspection  of  S  and  T. 

4.3.3  The  design  of  reduced  order  optimal  digital  LQG  controllers 

ft  is  our  expectation  that  the  reduced  order  controllers  should  maintain  as  much 
robustness  as  possible  when  the  performance  of  the  system  satisfies  the  design  require¬ 
ments.  Therfore,  a  two-step  robustness  recovery  procedure  will  be  applied  as  follows: 

( 1 )  Design  LQ-regulators  via  (4-26a)  for  the  full  order  system  For  all  frequen¬ 
cies  where  the  weights  W(e*‘^^)=H<l)(e^“^)B/p  are  much  larger  than  unity,  the  LQ-regu- 
lator  sensitivity,  complementary  sensitivity,  and  loop  transfer  function  have  the  following 
properties. 


(Ti  ((I+Kc<I>(e'‘^T)B)-ij  =  i/^. 

(4-27) 

(Ti  (Kc<I>(e'“T)B(HKcd)(e'‘"^B)"'|  =  I 

(4-28) 

(Ji  [KcdJ(ej‘‘’0B]  =  cTj  (W(e)"'')| 

(4-29) 

for  each  singular  value,  at  Therefore,  the  accuracy  may  be  improv'd  by  properly  in 
creasing  the  weight,  W(7.),  The  adjustible  parameters  are  H  and  p  lor  this  case;  if  H 
has  been  selected,  the  accuracy  may  be  improved  by  proper  reduction  of  the  parameter, 
p.  The  preliminary  design  of  the  LQ-regulator  just  invoices  the  selection  of  the  paramc 
ter,  p,  in  this  case. 

(2)  Design  a  sequence  of  Kalman  filters  via  (4-2 6b)  for  the  full  order  system 
Let  L=B,  and  p,  — ^  0,  to  approximate  the  function  in  step  1  to  whatever  robustness 

is  needed.  The  preliminary  design  of  the  Kalman  filter  involves  just  the  selection  of 
the  parameter,  p. 

(3)  Simulations  of  LQG  optimal  digital  control  for  various  reduced  order  control 
lers  for  the  .selected  parameter  pair  (p,p).  The  simulations  wall  be  used  for  verifying 
which  of  the  available  reduced  order  controllers  is  best. 


4.4  Numerical  Results 


It  is  assumed  that  the  plant  model  of  the  orbiting  shallow  spherical  shell  includes 
3  rigid  modes, 3  axisymmetric  modes,  I  meridional  mode,  6  combined  modes,  i.e.. 
26  dimensions  in  all.  The  controllers  may  be  divided  into  5  cases; 

Case  1  :  full  order  controllers 
Case  2  :  reduced  order  18-dim  controllers 
Case  3  ;  reduced  order  12-dim  controllers 
Case  4  :  reduced  order  8-dim  controllers 
Case  5  :  reduced  order  6-dim  controllers 

The  parameters  for  the  simulations  are  selected  as  follows: 

Sampling  time:  5  seconds 
System  noise  ;  Gs  =  10'^  (rad.) 

Observational  noise  :  Go  =  10"^  (meter) 

Physical  and  geometrical  parameters  of  the  shell: 
mass  =  10,000  kg. 
the  base  radius  of  shell  =  1  >0  meters 
the  height  of  shell  =  1  metei 
radius  of  curvature  for  shell  =  5000  meters 
wall  thickness  of  shell  =  0.01  meter 

Initial  conditions  for  all  simulatic  ns: 

\|/(0)=<|)(0)=e(0)=0.2  rad. 

4/(0)=<i)(0)=6(0)=0.02  rad./sjc 
qi(0)=q2(0)=q3(0)=...  =q«(0  =5  meters 
q7(0)=q8(0)=...  =qio(0)=0 
q,(0)=q2(0)=...  =q,o(0)=0 

(1)  The  design  of  the  desirable  LQ  regulators 

The  free  parameters  for  the  regulators  are  H  and  p.  If  we  let  H=C,  the  p  becomes 
the  only  free  parameter.  The  weighting  function  and  Eq.(27)  indicate  that  with  p  proper 
ly  reduced,  the  loop’s  errors  in  the  presence  of  commands  and  disturbances  can  be 
made  small.  But  p  cannot  be  too  small;  otherwise  the  errors  in  the  presence  of  ob.serva- 
tional  noise  will  increase  with  the  reduction  of  p.and  the  robustness  will  be  degraded, 
often  resulting  in  divergence  of  the  transient  responses  for  the  reduced  order  controllers. 
Fig.  4-3  certifies  this  point  and  shows  that  p=l  is  a  proper  value. 

(2)  The  design  of  the  desirable  Kalman  filler 

By  means  of  the  procedures  for  the  robustness  recovery,  as  we  know,  the  robust¬ 
ness  of  the  LQG  control  system  willbe  increased  when  p,  is  reduced;  but  the  p.  cannot 
be  too  small,  otherwise  the  accuracy  of  ihe  LOG  control  system  will  be  degraded.  The 
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design  of  Ihe  control  system  involves  a  compromise  between  robustness  and  accuracy. 
Fig. 4-4  certifies  tliis  point.  The  transient  responses  show  that  the  proper  p,  values  are 
O.l  or  0.01 

(3)  The  determination  of  the  reduced  order  LQG  controllers  in  the  presence 

of  unmodeied  dynamic  uncertainties. 

Simulations  have  been  conducted  for  the  LQG  control  (  plant  is  26  dimension;?) 
with  4  kinds  of  reduced  order  controllers(  18-dim,  12-dim,  8-dim. and  6-dim)  for 
the  orbiting  shallow  spherical  shell,  where  the  parameters,  p  and  p  are  set  as  1  and 
O.l,  respectively.  Fig. 4-5  shows  that  the  robustness  and  performance  of  the 
reduced  order  LQG  controllers  are  worse  than  that  of  the  full  order  controllers;  the 
12-dim  reduced  order  controller  is  sufficient  for  handling  the  unmodeled  dynamics 
of  the  shallow  spherical  shell  system;  the  6-dim  and  8-dim  reduced  order  controllers 
cannot  satisfy  the  requirement  of  the  optimal  control  of  the  shallow  shell  system;  the 
6-dim  controller,  which  is  based  only  on  the  rigid  modes,  will  result  in  severe  divergence 
for  the  transient  respon.ses  of  the  shallow  shell  system. 

4.5  Summary 

(1)  The  properties  of  robustness  recovery(sensitivit  recovery)  for  discrete-time 
systems  are  studied,  and  may  be  used  for  the  design  of  educed  order  optimal  LQG 
digital  controllers  forthe  shallow  spherical  shell  system.  The  design  of  the  control  system, 
in  fact,  is  a  compromise  between  robustness  and  performance  of  the  control  system 
after  the  .stability  of  the  control  system  is  satisfied  If  L=B  and  H=C,  the  robustness 
and  performance  of  the  control  system  only  depend  on  two  parameters:  p(sensitivity) 
and  p(robustnessi).  Therefore,  the  design  of  the  LQG  control  sy.stem  involves  tho  .selec¬ 
tion  of  the  proper  parameter  pair  (p,M.),  so  that  the  system  is  .stable  and  the  robustness 
and  performance  satisfy  the  design  requirements. 

(2)  If  the  robustness  recovery  is  used  for  increasing  the  robu.stness  of  the  control 
sy.stem,  in  general,  the  p  value  should  be  as  small  as  possible;  but  considering  the  per¬ 
formance  of  the  system,  the  p  cannot  be  selected  too  small,  otherwise  the  performance 
of  the  system  will  be  degraded.  In  general,  the  performance  of  the  reduced  order  control 
lers  is  worse  than  that  of  the  full  order  conn  oilers 

(3)  Simulations  have  certified  the  1 2-dim  reduced  order  controller  will  be 
sufficient  for  the  optimal  LQG  control  of  the  shallow  spherical  shell  system  in  the  pres 
ence  of  unmodeled  dynamics.  The  performance  of  the  8-dim  reduced  order  LQG  con 
trollers  for  the  shell  system  is  unacceptable.  The  6-dim  reduced  order  controllers  for 
the  shell  system  will  result  in  the  severe  divergence  of  the  transient  responses. 
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Fig. 4-3  The  Influence  of  p  in  the  Controller  on  the  Transient  Response  of 
LQG  Control  for  an  Orbiting  Shallow  Spherical  Shell  System 
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Fig. 4-5  The  Robustness  Comparison  of  Various  Reduced  Order  LQG  Controllers 
for  an  Orbiting  Shallow  Spherical  Shell  System  fp=l,  n=0.1) 
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5  The  Comparison  Between  the  Digital  Optimal  LQG  Controller  with  the 

Filtering  observer  and  Predicting  Observer  for  the  Orbiting  Flexibie  Shaliow 
Spherical  Shell  System 

As  we  know,  there  are  two  kinds  of  different  LQG  controllers  for  the  discrete-time 
system:  one  of  them  is  shown  Fig. 3-1,  in  which  the  state  being  fed  into  the  controller 
is  the  filtered  estimate  of  the  state  variable,  so  it  called  the  digital  LQG  controller  with 
filtering  observer;  another  is  shown  Fig. 3-2,  in  which  the  state  being  fed  into  the  control¬ 
ler  is  the  predicted  estimate  of  the  state  variable,  so  it  is  called  the  digital  LQG  controller 
with  the  predicting  observer. 

It  was  shown  in  Chapter  3  that  there  is  a  robustness  recovery  property  for  the 
digital  LQG  controller  with  filtering  observer,  but  there  is  no  corresponding  property 
for  the  digital  LQG  controller  with  the  predicting  observer.  The  proof  is  given  in  detail 
in  Chapter  3.  In  this  chapter,  the  comparison  between  the  transient  responses  of  the 
LQG  digital  control  with  filtering  observer  and  predicting  observer  will  be  made  by 
simulations  of  the  LQG  digital  control  for  the  orbiting  shallow  spherical  shell  sys- 
tem(Fig.4-l). 

It  is  assumed  that  the  plant  model  of  the  orbiting  shallow  spherical  shell  includes 
3  rigid  modes,  3  axisymmetric  modes,  I  meridional  mode,  6  combined  modes,  i.e., 
26  dimensions  ir  all.  The  controllers  may  be  divided  into  3  cases: 

Case  1  :  full  order  controllers  (3  rigid  modes  +  10  flexible  modes) 

Case  2  :  18-dim  reduced  order  controllers(3  rigid  modes  +  6  flexible  modes) 
Case  3  :  12-dim  reduced  order  controllers!  3  rigid  modes  +  3  flexible  modes) 

The  physical  and  geometrical  parameters  of  the  .<;hell: 
mass  =  10,000  kg. 
the  base  radius  of  shell  =100  meters 
the  height  of  shell  =  1  meter 
the  radius  of  curvature  for  shell  =  5000  meters 
wall  thickness  of  shell  =  0.01  meter 

The  parameters  for  the  simulations  are  selected  as  follows: 
sampling  time:  5  seconds 
system  noies:  a,  =  10"^  (rad.) 
observational  noise:  cTo  =  10"^  (meter) 

Initial  conditions  for  all  simulations: 
v)/(O)=<t)(O)=0(O)=O.2  rad. 

4/(0)=4)(0)=e(0)=0.02  rad. /sec. 
qi(0)=q2(0)=q3(0)=...=q6(0)  =  5  meters 
q7(0)=q8(0)=...=qio(0)=0.0 
qi(0)=q2(0)=...=q,o(0) 
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Based  on  the  development  in  the  Chapter  4,  the  relationship  between  the  robustness 
recovery  and  the  accuracy  of  the  control  system  may  be  listed  as  follows: 


Table  5-1  Robustness  recovery  and  system  accuracy 


Robustness  Recovery 

System  Accuracy 

Let  L  =  B,  r  ^  m 

Suppose  C<l>L  is  min.  phase 

where  <I>  =  (zI-A)"' 

Sensitivity  function 

S(z)  =  (l+T,(z))-' 

Sensitivity  weighting 

W(z)  =  Hd>B/p 

When  p,  0 

then  Ti(z)  Kc^B 

when  p,  0 

then  ai(S(z))— (Ti((  l+Ke<I)B)-*)=l/cTi(W) 

where 

Ti(z)=Gc(z)Gp(z)  :  the  transfer  function  of  the  input  loop-breaking  point  1 
Kc4)B  :  LQR  loop  transfer  ‘‘unction 

fn  the  design  of  LQG  r.ibust  control  systems,  the  free  parameters  are,  p,  the  control 
weighting  factor  in  performance  index  and,  p,,  the  measurement  noise  parameter.  As 
shown  in  Table  5-1,  the  reduction  of  p,  will  increase  the  system  robustness;  the  reduction 
of  p  will  reduce  the  ai(S(z)).  S(z)  is  the  sensitivity  function  of  the  sy.stem,  the  reduction 
of  CTi(S)  means  that  the  accuracy  of  the  control  system  will  be  improved. 


5.1  Comparison  of  transient  responses  between  the  full  order  LQG  digital 
controller  with  predicting  observer  and  filtering  observer 

In  order  to  study  and  compare  the  influence  of  ihe  .sensitivity  parameter,  p,  and 
robustness  parameter,  p,,  on  the  transient  response  of  the  digital  full  order  LQG  control¬ 
ler  with  predicting  observer  and  filtering  observer,  the  parameters  p,  and  p.  are  varied 
in  the  combination  shown  in  Table  5-2  and  Table  5-3  for  the  simulations: 


Table  5-2  The  parameter  pair  (p,p,)  in  Fig.  5-1 


Case 

Predicting  observer 

Filtering  observer 

1-1 

1 

1 

Fig.5-1-1 

Fig.5-1-2 

1-2 

1 

0.1 

Fig.5-1-3 

Fig. 5-1-4 

1-3 

1 

0.01 

Fig.5-1-5 

Fig.5-1-6 

Table  5-3  The  parameter  pair  (p.p.)  in  Fig.  5-2 


Case 

M- 

Predicting  observer 

Filtering  observer 

1-4 

0.1 

1 

Fig.  5-2-1 

Fig.  5-2-2 

1-5 

0.1 

0.1 

Fig.5-2-3 

Fig.  5-2-4 

1-6 

0.1 

0.01 

Fig.5-2-5 

Fig.  5-2-6 

The  following  two  points  are  shown  in  Fig. 5-1  and  Fig.5-2: 

1  Fig. 5-1  shows  that  the  transient  respon.se  of  the  digital  full  order  LQG  control¬ 
ler  with  filtering  observer  for  the  orbiting  shallow  spherical  shell  system  is  more  sensitive 
to  the  reduction  of  the  parameter,  p,,  than  that  of  the  digital  full  order  LQG  controller 
with  the  predicting  observer.  The  transient  performance  of  the  system  with  filtering  ob¬ 
server  is  better  than  that  of  the  system  with  predicting  observer  for  the  same  parameter 
pair  (p,^t). 

2  Fig.5-2  shows  that  the  decay  of  the  transient  response  of  the  system  will  in¬ 
crease  strongly  with  the  reduction  of  the  sensitivity  parameter,  p;  i.e.,  the  time  of  the 
transient  process  of  system  will  be  shortened  when  the  value  of  p  is  reduced.  In  this 
case  the  phenomenon  of  the  data  saturatim  will  ilso  appear  earlier  in  the  response, 
but  the  accuracy  of  the  system  with  filtering  ohscr\cr  is  still  better  than  the  system  with 
predicting  observer. 


5.2  Comparison  of  transient  responses  between  the  18-dim  reduced  order  LQG 
digital  controller  with  predicting  observer  and  filtering  observer 

In  order  to  study  and  compare  the  influence  of  the  sensitivity  parameter,  p.  and 
robustness  parameter,  p.,  on  the  transient  response  of  the  digital  18-dim  reduced  order 
LQG  controller  with  predicting  observer  and  filtering  observer,  the  parameters  p,  and 
p,  are  varied  in  the  combinations  shown  in  the  Tables  5-4  and  5-5  for  the  simulations: 
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Table  5-4  The  parameter  pair  (p,|x)  in  Fig.  5-3 


Case 

inQiiim 

Predicting  observer 

Filtering  observer 

2-1 

1 

1 

Fig.5-3-1 

Fig.5-3-2 

2-2 

1 

0.1 

Fig.5-3-3 

Fig. 5-3-4 

2-3 

1 

0.01 

Fig.5-3-5 

Fig.5-3-6 

Table  5-5 

The  parameter  pair  (p,p)  in  Fig.  5-4 

Case 

P 

M- 

Predicting  observer 

Filtering  observer 

2-4 

0.1 

1 

Fig.5-4-1 

Fig.  5-4-2 

2-5 

0.1 

0.1 

Fig. 5-4-3 

Fig.  5-4-4 

2-6 

0.1 

0.01 

Fig.  5-4-5 

Fig.5-4-6 

The  following  points  are  shown  in  Fig. 5-3,  and  Fig. 5-4: 

1  Fig. 5-3  shows  that  the  18-dim.  reduced  order  LQG  digital  controller  for  the 
both  types  of  observer  are  stable  for  the  parameters  p=l  and  0. 1,  0.01;  The  tran¬ 
sient  response  of  the  system  will  be  improved  when  the  value  of  the  parameter,  p.  is 
decreased  from  I  to  0.01.  The  performance  of  the  system  with  the  filtering  observer 
is  better  than  that  of  the  system  with  the  predicting  observer  when  the  value  of  the  param¬ 
eter,  p,  is  reduced. 

2  Fig. 5-4  shows  that  the  transient  response  of  the  system  with  the  filtering 
observer  for  the  orbiting  shallow  -pherical  shell  system  is  still  stable  when  the  values 
of  the  parameter,  p,  is  reduced,  but  the  transient  response  of  the  system  with  the  predict¬ 
ing  observer  is  not  stable  for  the  parameters  p=0. 1,  p=0.01.  This  result  is  due  to  the 
fact  that  the  sensitivity  of  the  system  will  be  increased  when  the  value  of  the  parameters, 
P,  is  decreased  from  1  to  0. 1 .  This  indicates  that  the  robustness  of  the  system  is  reduced 
with  the  reduction  of  the  parameter,  p.  Therefore,  the  transient  response  of  the  system 
with  predicting  observer  will  result  in  divergence  since  the  robustness  of  system  is  not 
sufficient  to  cover  the  error  of  the  un modelled  dynamics;  but  the  transient  response 
of  the  system  with  the  filtering  observer  is  still  stable  due  to  the  robustness  recovery 
property  of  the  system  with  filtering  observer. 


5.3  Comparison  of  transient  responses  between  the  12-dim  reduced  order  LQG 
digital  controller  with  predicting  observer  and  filtering  observer 

In  order  to  study  and  compare  the  influence  of  the  sensitivity  parameter,  p,  and 
the  robustness  parameter,  p.  on  the  transient  response  of  the  digital  12-dim.  reduced 
order  LQG  controller  with  the  predicting  observer  and  filtering  observer,  the  paramcici  s. 
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p,  and  VJL  are  combined  indicated  in  Tables  5-6  and  5-7  for  the  simulations: 


Table  5-6  The  parameter  pair  (p.p.)  in  Fig.  5-5 


Case 

M- 

Predicting  observer 

Filtering  observer 

3-1 

1 

1 

Fig.5-5-1 

Fig.  5-5-2 

3-2 

1 

0.1 

Fig. 5-5-3 

Fig. 5-5-4 

3-3 

1 

0.01 

Fig.5-5-5 

Fig.  5-5-6 

Table  5- 

-7  The  parameter  pair  (p,p)  in  Fig.  5-6 

Case 

M- 

Predicting  observer 

Filtering  observer 

3-4 

0.1 

1 

Fig.  5- 6-1 

Fig.  5-6-2 

3-5 

0.1 

0.1 

Fig.5-6-3 

Fig.  5-6-4 

3-6 

0.1 

0.01 

Fig.  5-6-5 

Fig. 5-6-6 

The  following  points  are  shown  in  the  Fig  5-5  and  Fig.  5-6: 

1  Fig.  5-5  shows  that  the  12-dim.  reduced  order  LQG  digital  controller  for 
both  observers  are  stable  for  p  =  I,  and  |ji=l,0. 1,0.01.  The  transient  response  of  the 
system  will  be  improved  when  the  value  of  p,  decreases  from  1  to  0.0 1.  The  performance 
of  the  sy.stem  with  the  filtering  observer  is  better  than  that  of  the  system  with  the  predicting 
observer  when  the  value  of  p-  is  reduced. 

2  Fig.  5-5  and  Fig.  5-6  show  that  transient  response  of  the  systems  with  predict¬ 
ing  observer  and  filtering  observer  are  not  stable,  because  the  system  robustness  is  re 
duced  when  the  value  of  p  is  reduced  from  1  to  ').01,  and  the  error  of  unmodelled 
dynamics  for  the  12-dim.  reduced  order  controller  is  greater  than  that  for  the  18-dim. 
reduced  order  controller,  the  robustness  of  these  sy.stems  is  not  sufficient  to  overcome 
the  unmodelled  dynamics.  Therefore,  the  divergence  of  the  system’s  transient  response 
results  for  all  parameter  pairs  (p,p)  in  Fig.  5-6. 

What  we  should  point  out  is  that  the  divergence  of  the  system  with  the  filtering 
observer  is  slower  than  that  of  the  system  with  the  predicting  observer. 
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Fig.5-2  The  Comparison  Between  the  Modal  Amplitude  Response  of  the  Full  Order 
LQ(3  Controller  with  Predicting  Observer  and  Filtering  Observer  (p=0.1) 
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Fig.5-3  The  Comparison  Between  the  .Modal  Amplitude  Response  of 
the  18-Dim.  Reduced  Order  LQG  Controller  with  Predicting 
Observer  and  Filtering  Observer  (p=l) 
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Fig.  5-4  The  Comparison  Between  the  Modal  Amplitude  Response  of 
the  18-Dlm.  Reduced  Order  LQG  Controller  with  Predicting 
Observer  and  Filtering  Observer  (p=0.1) 
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Fig.  5-5  The  Comparison  Between  the  Modal  Amplitude  Response  of 
the  1 2-Dim  Reduced  Order  LQG  Controller  with  Predicting 
Observer  and  Filtering  Observer  (p-l) 
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6  General  Conclusions  and  Suggested  Future  Direction 
6.1  General  conclusions 

In  order  to  develop  the  analysis  and  design  methods  for  robustness  control  of  large 
space  structure  sampled  data  stochastic  systems  with  a  specific  application  to  the  orbiting 
flexible  shallow  spherical  system,  the  following  problems  have  been  studied  and  the 
conclusions  we  have  obtained  are  as  follows: 

(1)  The  theory  of  multi-input,  multi-output(MIMO)  transfer  function  matrices 
in  the  z  transformation  may  be  used  for  the  analysis  and  design  in  the  frequency  domain 
for  the  discrete-time  system.  The  frequency  response  of  the  transfer  function  matrix 
for  the  discrete-time  system  can  be  obtained  by  means  of  the  frequency  response  of 
the  transfer  function  matrix  for  the  equivalent  continous-time  system. 

(2)  The  robustness  criteria  in  the  frequency  domain  for  discrete-time  systems 
have  been  developed.  The  stability  conditions  for  MIMO  discrete-time  feedback  system 
with  additive  alteration  are  (2-9)  or  (2-10);  The  stability  conditions  for  MIMO  dis¬ 
crete-time  feedback  with  multiplicative  alteration  are  (2-14)  or  (2-16). 

(3)  One  of  the  most  important  breakthroughs  in  multi-input,  multi-output  feed¬ 
back  system  theory  for  the  last  decade  is  the  development  of  the  loop  transfer  recovery 
methodology  for  the  continous-time  linear  quadractic  Gaussian  problem,  which  is  called 
LQG/LTR.  The  LQG/LTR  technique  has  been  extended  from  the  continous-time  sys¬ 
tem  to  the  discrete-time  system  in  this  report.  It  is  proven  that  the  robusmess  (sensitivity) 
recovery  property  is  also  valid  for  the  LQG  digital  controller  with  the  filtering  observer, 
but  it  is  not  valid  for  the  LQG  digital  controller  with  the  prediaing  observer. 

(4)  As  an  application  of  the  LQG/LTR  technique  for  discrete-time  systems  to 
large  space  struaural  systems,  the  LQG/LTR  technique  of  the  discrete-time  systems  has 
been  used  to  design  the  reduced  order  optimal  digital  LQG  controller  for  the  orbiting 
flexible  shallow  spherical  shell  system.  The  research  results  indicate  that  the  12-dim 
reduced  order  controller  will  be  sufficient  for  the  optimal  LQG  control  of  the  shallow 
spherical  shell  system  in  the  presence  of  unmodelled  dynarrics. 

(5)  The  comparisons  between  the  digital  optimal  LQG  controller  with  the  filtering 
observer  and  predirting  observer  for  the  orbiting  flexible  shallow  spherical  shell  system 
have  been  made.  The  robustness  recovery  property  for  the  digital  LQG  controller  with 
filtering  observer  has  been  certified  by  the  simulations.  The  :>imulations  indicate  that 
the  transient  response  of  the  digital  LQG  control  system  with  the  filtering  observer  or 
with  the  predicting  observer,  in  general,  depends  on  the  robustness  parameter,  |x,  and 
sensitivity  parameter,  p.  The  best  comb.'nation  of  the  parameters  p  and  n,  will  depend 
on  the  compromise  between  the  accuracy  and  robustness.  Censidering  the  robustness 
recovery  property,  the  s>stem  performance  of  the  robust  control  system  with  the  filtering 
observer  will  be  better  than  that  of  the  LQG  robust  control  system  with  the  predicting 
observer. 


6.2  The  suggested  Future  directions 

(1)  Loop  transfer  recovery  for  nonminimum  phase  plants  for  discrete-time 
systems 

The  requirement  of  minimum  phase  plant  (i.c.,  the  transfer  function  of  the  plant 
has  no  (finite)  zero  outside  the  unit  circle)  for  the  recovery  procedure  is  critical.  Since 
there  are  some  plants  which  are  minimum  phase  systems  in  practice,  it  would  be  desir¬ 
able  to  have  a  methodology  for  incorporating  limitations  due  to  non-minimum  phase 
zeros  into  the  LTR  procedure.  It  is  especially  more  desirable  for  discrete-time  systems, 
since  the  standard  sampling  process  is  known  to  introduce  zeros,  some  of  which  some¬ 
time  lie  outside  the  unit  circle. 

(2)  Synthesis  and  design  of  reduced  order  LQG/LTR  optimal  digital 
controllers  using  constrained  optimization  techniques 

It  is  well  known  that  the  basic  requirements  of  a  feedback  system  are  better  per- 
formance(small  error  in  the  presence  of  disturbances  and  reference  input)  and  robust- 
ness(stability  and  performance  maintained  in  the  presence  of  model  uncertainties).  In 
fact,  the  two  parts  of  these  basic  requirements  are  in  conflict  with  each  other.  As  far 
as  the  LQG/LTR  method  is  concerned,  the  conflict  is  reflected  in  the  selections  of  the 
robustness  parameter  ai  d  sensitivity  parameierjSj.  Therefore,  we  may  convert  the  prob¬ 
lem  of  reduced  order  LQG/LTR  controller  design  into  the  constrained  optimization 
problem.  This  procedure  minimizes  a  linear  quadratic  Gaussian(LQG)  type  cost  func¬ 
tion  while  trying  to  satisfy  a  set  of  constraints  on  the  responses  and  stability  margins. 
Although  a  linear  LQG  cost  function  was  minimized  by  updating  the  free  parameters 
of  the  control  law,  while  satisfying  a  ^et  of  constraints  on  the  design  loads,  responses, 
and  stability  margin[18],  our  attention  will  Itc  focused  on  the  design  and  synthesis  of 
the  reduced  order  LQG/LTR  optimal  digital  controller  for  discrete-time  systems,  using 
only  a  small  number  of  design  parameters  spccifif  ally  associated  with  robustness  and 
sensitivity.  As  an  application,  this  method  will  Ik  applied  to  the  design  of  digital  reduced 
order  LQG/LTR  controllers  for  the  orbiting  shallow  spherical  shell  system. 

(3)  The  synthesis  and  design  of  the  robust  digital  optimal  output  feedback 
reduced  order  controllers  using  constrained  optimization 

Since  the  mathematical  system  model  is  inherently  of  high  order  for  large  space 
structural  systems  and  because  of  the  practical  possibility  of  on-board  computational 
implementation,  it  is  desirable  to  have  methods  avaible  for  the  design  of  low-order 
controllers  for  high-order  plants.  Such  methods  can  be  broadly  divided  into  two  classes: 
(a)  direct  method:  in  which  the  parameters  defining  a  low-order  controller  are  com 
puted  by  some  optimization  or  other  procedure;  (b)  indirect  method:  in  which  a  high- 
order  controller  is  first  found  and  then  a  procedure  used  to  simplify,  or  a  low-order 
plant  first  is  found  by  some  criterion,  and  then  a  low-order  controller  based  on  the 
simplified  low-order  plant  is  designed.  In  general,  the  direct  method  is  better  than  the 
indirect  method  in  meeting  the  requirement  of  the  designer.  The  design  method  for  the 
robust  digital  optimal  output  feedback  reduced  order  controller  using  constrained  opti 
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mization  is  just  the  direct  design  method  for  the  low-order  controller.  Therefore,  it  is 
very  useful  to  study  and  develop  the  design  method  for  the  robust  digital  optimal  output 
feedback  reduced  order  controller  using  constrained  optimization. 
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